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Abstract 



The form factor ratios ry = V(0)/Ai(0), r 2 = A 2 (0)/A 1 (0) and r 3 = A 3 (0)/Ai(0) in 
the decay D + — > K l + ve, K — ► K~ir + have been measured using data from charm 
hadroproduction experiment E791 at Fermilab. From 3034 (595) signal (background) 
events in the muon channel, we obtain ry = 1.84 ± 0.11 ± 0.09, r 2 = 0.75 ± 0.08 ± 0.09 
and, as a first measurement of r 3 , we find 0.04 ± 0.33 ± 0.29. The values of the form 
factor ratios ry and r 2 measured for the muon channel are combined with the values of 
ry and r 2 that we have measured in the electron channel. The combined E791 results for 
the muon and electron channels are ry = 1.87 ± 0.08 ± 0.07 and r 2 = 0.73 ± 0.06 ± 0.08. 
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The weak decays of hadrons containing heavy quarks are substantially influ- 
enced by strong interaction effects. Semileptonic charm decays such as D + — > 
U£ are an especially clean way to study these effects because the leptonic 
and hadronic currents completely factorize in the decay amplitude. All informa- 
tion about the strong interactions can be parametrized by a few form factors. Also, 
according to Heavy Quark Effective Theory, the values of form factors for some 
semileptonic charm decays can be related to those governing certain b-quark de- 
cays. In particular, the form factors studied here can be related to those for the 
rare 5-meson decays B — > K*e + e~ and B — > K*^f |l],|2| which provide windows 
for physics beyond the Standard Model. 

With a vector meson in the final state, there are four form factors, V(q 2 ), 
Ai(q 2 ), A 2 (q 2 ) and A 3 (q 2 ), which are functions of the Lorentz- invariant momentum 
transfer squared ||. The differential decay rate for D + — > K* C p^v^ with K*° — » 
K~7i + is a quadratic homogeneous function of the four form factors. Unfortunately, 
the limited size of current data samples precludes precise measurement of the (in- 
dependence of the form factors; we thus assume the dependence to be given by 
the nearest-pole dominance model: F(q 2 ) = F(0)/(1 — q 2 /m 2 ole ) where m po i e = 
my = 2.1 GeV/c 2 for the vector form factor V, and m po i e = tua = 2.5 GeV/c 2 
for the three axial- vector form factors [|J. The third form factor A 3 (q 2 ), which is 
unobservable in the limit of vanishing lepton mass, probes the spin-0 component of 
the off-shell W. Additional spin-flip amplitudes, suppressed by an overall factor of 
m 2 /q 2 when compared with spin no-flip amplitudes, contribute to the differential 
decay rate. Because Ai(q 2 ) appears among the coefficients of every term in the 
differential decay rate, it is customary to factor out Ai(0) and to measure the 
ratios r v = V(0)/Ax(0), r 2 = A 2 {0)/A 1 (0) and r 3 = A 3 {a)/A l (0). The values of 
these ratios can be extracted without any assumption about the total decay rate 
or the weak mixing matrix element V cs . 

We report new measurements of the form factor ratios for the muon chan- 
nel and combine them with slightly revised values of our previously published 
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measurements of ry and r2 f| for the electron channel. This is the first set of 
measurements in both muon and electron channels from a single experiment. We 
also report the first measurement of r 3 = A 3 (0)/Ai(0), which is unobservable in 
the limit of vanishing charged lepton mass. 

E791 is a fixed-target charm hadroproduction experiment 0. Charm particles 
were produced in the collisions of a 500 GeV/c 7r _ beam with five thin targets, 
one platinum and four diamond. About 2 x 10 10 events were recorded during the 
1991-1992 Fermilab fixed-target run. The tracking system consisted of 23 planes 
of silicon microstrip detectors, 45 planes of drift and proportional wire chambers, 
and two large-aperture dipole magnets. Hadron identification is based on the in- 
formation from two multicell Cerenkov counters that provided good discrimination 
between kaons and pions in the momentum range 6 — 36GeV/c. In this momentum 
range, the probabability of misidentifying a pion as a kaon depends on momentum 
but does not exceed 5%. We identified muon candidates using a single plane of 
scintillator strips, oriented horizontally, located behind an equivalent of 2.4 meters 
of iron (comprising the calorimeters and one meter of bulk steel shielding). The an- 
gular acceptance of the scintillator plane was Rdz62mrad x ±48mrad (horizontally 
and vertically, respectively), which is somewhat smaller than that of the rest of the 
spectrometer for tracks which go through both magnets (~ ±100mrad x ±64mrad). 
The vertical position of a hit was determined from the strip's vertical position, and 
the horizontal position of a hit from timing information. 

The event selection criteria used for this analysis are the same as for the 
electronic-mode form factor analysis || , except for those related to lepton identifi- 
cation. Events are selected if they contain an acceptable decay vertex determined 
by the intersection point of three tracks that have been identified as a muon, a 
kaon, and a pion. The longitudinal separation between this candidate decay vertex 
and the reconstructed production vertex is required to be at least 15 times the esti- 
mated error on the separation. The two hadrons must have opposite charge. If the 
kaon and the muon have opposite charge, the event is assigned to the "right-sign" 
sample; if they have the same charge, the event is assigned to the "wrong-sign" 
sample used to model the background. 

To reduce the contamination from hadron decays in flight, only muon can- 
didates with momenta larger than 8 GeV/c are retained. With this momentum 
restriction, the efficiency of muon tagging was about 85%, and the probability for 
a hadron to be identified as a muon was about 3%. 

To exclude feedthrough from D + — > K~7i + 7r + , we exclude events in which the 
invariant mass of the three charged particles (with the muon candidate interpreted 
as a pion) is consistent with the D + mass. For our final selection criteria, we use 
a binary-decision-tree algorithm (CART [^]), which finds linear combinations of 
parameters that have the highest discrimination power between signal and back- 
ground. Using this algorithm, we found a linear combination of four discrimination 
variables || : (a) separation significance of the candidate decay vertex from target 
material; (b) distance of closest approach of the candidate D momentum vector 
to the primary vertex, taking into account the maximum kinematically-allowed 
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miss distance due to the unobserved neutrino; (c) product over candidate D decay 
tracks of the distance of closest approach of the track to the secondary vertex, 
divided by the distance of closest approach to the primary vertex, where each dis- 
tance is measured in units of measurement errors; and (d) significance of separation 
between the production and decay vertices. This final selection criterion reduced 
the number of wrong-sign events by 50%, and the number of right-sign events by 
25%. Although this does not affect our sensitivity substantially, it does reduce 
systematic uncertainties associated with the background subtraction. 

The minimum parent mass M m i n is defined as the invariant mass of Knfiv when 
the neutrino momentum component along the D + direction of flight is ignored. The 
distribution of M m ; n should have a Jacobian peak at the D + mass, and we observe 
such a peak in our data (Fig. 1). We retain events with M min in the range 1.6 
to 2.0 GeV/c 2 as indicated by the arrows in the figure. The distribution of Ktt 
invariant mass for the retained events is shown in the top right of Fig. 1 for both 
right-sign and wrong-sign samples. Candidates with 0.85 < Mk-k < 0.94GeV/c 2 
were retained, yielding final data samples of 3629 right-sign and 595 wrong-sign 
events. 

The hadroproduction of charm, the differential decay rate, and the detector 
response were simulated with a Monte Carlo event generator. A sample of events 
was generated according to the differential decay rate (Eq. 22 in Ref. Pf), with the 
form factor ratios ry = 2.00, r2 = 0.82, and r% = 0.00. The same selection criteria 
were applied to the Monte Carlo events as to real data. Out of 25 million generated 
events, 95579 decays passed all cuts. Figure 1 (bottom) shows the distribution of 
Mktt from real data after background subtraction ("right-sign" minus "wrong- 
sign") overlaid with the corresponding Monte Carlo distribution after all cuts are 
applied. The agreement between the two distributions suggests that wrong-sign 
events correctly account for the size of the background. 

The differential decay rate |5| is expressed in terms of four independent kine- 
matic variables: the square of the momentum transfer (q 2 ), the polar angle 9y in 

the K*° rest frame between the 7r + and D + , the polar angle 9g. in the W + rest 
frame between the and D + , and the azimuthal angle x i n the D + rest frame 
between the K and W + decay planes. The definition we use for the polar angle 
6 1 is related to the definition used in Ref. |§ by 9e — > it — 9e. 

Semileptonic decays cannot be fully reconstructed due to the undetected neu- 
trino. With the available information about the D + direction of flight and the 
charged daughter particle momenta, the neutrino momentum (and all the decay's 
kinematic variables) can be determined up to a two-fold ambiguity if the parent 
mass is constrained. Monte Carlo studies show that the differential decay rate is 
more accurately determined if it is calculated with the solution corresponding to 
the lower laboratory-frame neutrino momentum. 

To extract the form factor ratios the distribution of the data points in the four- 
dimensional kinematic variable space is fit to the full expression for the differential 
decay rate. We use the same unbinned maximum-likelihood fitting technique as in 
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our D + —> K e + z/ e form factor analysis JS|. The likelihood function is computed 
from the density of weighted Monte Carlo events (described above) near each 
data event in the four-dimensional space of kinematic variables ||. To include 
background in the fit, a similar likelihood function based on the density of wrong- 
sign events around each right-sign event is used. With this method the fitted 
results are subject to small systematic biases which originate from two sources: 
(a) approximate normalization of the likelihood function; (b) nonlinearity of the 
decay rate within the volume centered on the data point. These systematic biases 
of the fitted parameters were determined from Monte Carlo studies, and are 5ry = 
+0.09 ± 0.02, 5r 2 = -0.08 ± 0.01 and 5r 3 = -0.11 ± 0.06. After correction for 
these biases by subtracting these 5r's from the measured values, the final form 
factor ratios and their statistical errors are ry = 1.84 ± 0.11, r 2 = 0.75 ± 0.08 
and r 3 = 0.04 ± 0.33. The sensitivity of the fit to r 3 is low because this form 
factor ratio only contributes to spin-flip amplitudes. The correlation coefficient 
between the form factors to which we are most sensitive, ry and r 2 , is —0.090. 
The correlation coefficient between r 3 and r 2 is —0.211, and between r 3 and ry 
is —0.087. Figure 2 compares the data and Monte Carlo distributions in various 
regions of the four- dimensional phase space. 

We checked for any potential bias in these results due to our choice of neu- 
trino momentum by employing a secondary fitting technique. Again, Monte Carlo 
is used to account for detector acceptance and smearing. However, both solu- 
tions for the neutrino momentum are now used in the fit. We divide the four- 
dimensional kinematic variable space into 240 separate volumes and determine the 
number of data entries in each volume, where each event has two entries - one 
for each neutrino-momentum solution. We use Monte Carlo events to determine 
the probabilities that an event generated in a particular phase space volume will 
be observed in each of the other volumes when the wrong neutrino-momentum 
choice is used. This feedthrough probability matrix and the observed number of 
data events determine the fraction of data events that correspond to the correct 
neutrino-momentum solution in each volume of kinematic variable space. Each 
fraction is then used in a binned maximum likelihood fit. Background is modeled 
with wrong-sign candidates as in the primary method. The form factor ratios 
and statistical errors measured with this secondary method are ry = 1.90 ± 0.11, 
r 2 = 0.72 ± 0.08, and r 3 = —0.25 ± 0.34. This method for extracting the form 
factor ratios uses the same D + — > K [i^v^ candidates as the previous method, 
but uses additional neutrino-momentum solutions. This is true for both the data 
and for the Monte Carlo sample used in the likelihood function calculation, so the 
results of this fit could differ from those of the previous fit. 

The values of the form factor ratios obtained with the two methods agree well, 
providing further assurance that selecting the lower neutrino momentum solution 
in the primary method and correcting for the systematic bias gives the correct 
result. However, the systematic uncertainties for the primary method (see below) 
were found to be significantly smaller, mainly because the unbinned maximum- 
likelihood method is more stable against changes in the size of the phase space 
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volume. Therefore, the primary method was chosen for quoting final results. 

We classify systematic uncertainties into three categories: (a) Monte Carlo 
simulation of detector effects and production mechanism; (b) fitting technique; (c) 
background subtraction. The estimated contributions of each are given in Table |. 
The main contributions to category (a) are due to muon identification and data 
selection criteria. The contributions to category (b) are related to the limited size 
of the Monte Carlo sample and to corrections for systematic bias. 

The measurements of the form factor ratios for D + — > K*° fi + u^ presented here 

and for the similar decay channel D + —> K e + v e || follow the same analysis 
procedure except for the charged lepton identification. Both results are listed in 
Table |T|. The consistency within errors of the results measured in the electron 
and muon channels supports the assumption that strong interaction effects, incor- 
porated in the values of form factor ratios, do not depend on the particular W + 
leptonic decay. Based on this assumption, we combine the results measured for 
the electronic and muonic decay modes. The averaged values of the form factor 
ratios are r v = 1.87 ± 0.08 ± 0.07 and r 2 = 0.73 ± 0.06 ± 0.08. The statistical and 
systematic uncertainties of the average results were determined using the general 
procedure described in Ref. M (Eqns. 3.40 and 3.40'). Some of the systematic 
errors for the two samples have positive correlation coefficients, and some nega- 
tive. The combination of all systematic errors is ultimately close to that which one 
would obtain assuming all the errors are uncorrelated. The third form factor ratio 
r 3 was not measured in the electronic mode. 

Table |T| compares the values of the form factor ratios ry and r 2 measured 
by E791 in the electron, muon and combined modes with previous experimental 
results. The size of the data sample and the decay channel are listed for each case. 
All experimental results are consistent within errors. The comparison between the 
E791 combined values of the form factor ratios ry and r 2 and previous experimental 
results is also shown in Fig. 3 (top). Table [TTTI and Fig. 3 (bottom) compare the final 
E791 result with published theoretical predictions. The spread in the theoretical 
results is significantly larger than the E791 experimental errors. 

To summarize, we have measured the values of the form factor ratios in the 
decay channel D+ -> 7T*V + ^ to be r v = 1.84 ± 0.11 ± 0.09, r 2 = 0.75 ± 0.08 ± 
0.09 and r 3 = 0.04 ± 0.33 ± 0.29. The data sample has about 3000 events after 
subtracting 595 background events. Combining these results for ry and r 2 with 
those measured by E791 for the decay channel D + — > K e + i> e gives ry = 1.87 ± 
0.08 ± 0.07 and r 2 = 0.73 ± 0.06 ± 0.08. 
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Academy of Sciences and Humanities, the U.S. Departament of Energy, the U.S.- 
Israel Binational Science Foundation, and the U.S. National Science Foundation. 
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TABLES 



TABLE I. The main contributions to uncertainties on the form factor ratios. 



Source 


0Y 2 


a ry 


o> 3 


Simulation of detector effects: 


0.06 


0.08 


0.16 


Hadron identification 


0.01 


0.01 


0.02 


Muon identification 


0.04 


0.06 


0.10 


Production mechanism 


0.01 


0.01 


0.02 




01 


02 


08 


Cut selection 


0.03 


0.04 


0.09 


Fitting technique: 


0.02 


0.03 


0.22 


MC volume size 


0.02 


0.02 


0.12 


Number of MC points 


0.01 


0.01 


0.18 


Bias 


0.01 


0.02 


0.06 


Background: 


0.06 


0.04 


0.08 


No. of background events 


0.04 


0.02 


0.06 


Background shape 


0.04 


0.04 


0.06 


Total 


0.09 


0.09 


0.29 



TABLE II. Comparison of E791 results with previous experimental results. The 
E791 electron result for ry is 0.06 higher than the value reported in Ref. [5] because we 
have corrected for inaccuracies in the earlier modeling of the D + transverse momentum. 



Exp. 



Events 



r v = V(0)/A 1 (0) 



r 2 = A 2 (P)/A 1 (p) 



E791 
E791 
E791 
E687 II 



E653 111 



E691 (T§ 



6000 (e + n) 
3000 (n) 
3000 (e) 
900 (n) 
300 (/J) 
200 (e) 



1.87 ±0.08 ±0.07 
1.84 ±0.11 ±0.09 
1.90 ±0.11 ±0.09 
1.74 ±0.27 ±0.28 



2.00 



+0.34 
-0.32 



±0.16 



2.0 ±0.6 ±0.3 



0.73 ± 0.06 ± 0.08 
0.75 ± 0.08 ± 0.09 
0.71 ± 0.08 ± 0.09 
0.78 ±0.18 ±0.10 
0.82±g;|§ ±0.11 
0.0 ±0.5 ±0.2 
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TABLE III. Comparison of E791 results with theoretical predictions for the form 
factor ratios ry and ri- 



Group 


r v 


T2 


E791 (e and fi) 


1.87 ±0.11 


0.73 ±0.10 



ISGW2 |3J 2.0 1.3 

WSB 1.4 1.3 

KS @ 1.0 1.0 

AW/GS U|]| 2.0 0.8 

Stech 1.55 1.06 

BKS pi 1.99 ±0.22 ±0.33 0.70 ± 0.16 ± 0.17 

LMMSgl 1.6 ±0.2 0.4 ±0.4 



ELC ||g 1.3 ±0.2 0.6 ±0.3 



APE [E1J1 1.6 ±0.3 0.7 ±0.4 



UKQCD Q 1.4±gj 0.9 ±0.2 

BBD H| 2.2 ±0.2 1.2 ±0.2 

LANL [Oil 1.78 ±0.07 0.68 ±0.11 
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FIGURES 
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1.05 
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FIG. 1. Distributions of minimum parent mass M m i n and Ktt invariant mass for 
D + — > .fT candidate events. Right-sign (RS) and wrong-sign (WS) samples are 

denned in the text. Top left: M m - m for events with if7r mass in the range 0.85 to 
0.94 GeV/c 2 . Top right: Ktt invariant mass for events with M m i n in the range 1.6 to 
2.0 GeV/c 2 . Bottom: background-subtracted (RS-WS) Kit mass distribution (crosses) 
compared to Monte Carlo prediction (dashed histogram) for events with M m i n in the 
range 1.6 to 2.0 GeV/c 2 . All candidates pass all the other final selection cuts. The 
arrows indicate the range of the final sample. 
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FIG. 2. Comparison of single- variable distributions of background-subtracted data 
(crosses) with Monte Carlo predictions (dashed histograms) using best-fit values for the 
form factor ratios. 
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FIG. 3. Top: Comparison of experimental measurements of form factor ratios ry 
and r2 for D + — > K t^vg in the muon ( / u), electron (e) and combined (/x + e) channels. 
The smaller error bars indicate the statistical errors and the larger ones indicate the 
statistical and sytematic errors added in quadrature. Bottom: Comparison of theoretical 
predictions with the E791 [p, + e) result. 
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